Non-equilibrium delocalization-localization transition of photons in circuit QED 
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We show that photons in two tunnel-coupled microwave resonators each containing a single superconduct- 
ing qubit undergo a sharp non-equilibrium delocalization-localization (self-trapping) transition due to strong 
photon-qubit coupling. We find, that self-trapping of photons in one of the resonators (spatial localization) 
forces the qubit in the opposite resonator to remain in its initial state (energetic localization). This allows for an 
easy experimental observation of the transition by local read-out of the qubit state. Dissipation of photons and 
decoherence of the qubit favor the self-trapped regime. 



In circuit quantum electrodynamics (QED), superconduct- 
ing qubits are coupled with microwave photons in a transmis- 
sion line resonator reaching extremely strong light-matter in- 
teractions within an integrated circuit jit]. Device integration, 
high tunability and individual addressability of each resonator 
make wide parameter regimes easily accessible. Circuit QED 
thus constitutes one of the most promising solid-state architec- 
tures for quantum information processing and offers the pos- 
sibility to study fundamental questions of interacting quan- 
tum systems 121]. The experimental focus has been on the 
design of the coupling between a single cavity and a single 
qubit and subsequent work demonstrated a great level of con- 
trol of single-cavity systems Today, a key challenge 
for scalability and further progress in the field is the under- 
standing of small coupled systems, i.e., effective qubit-qubit 
and photon-photon interactions and their interplay with dissi- 
pation IgjjTili. 

Recent theoretical interest in hybrid light-matter systems 
has been on the superfluid-Mott insulator transition of polari- 
tons in a coupled-cavity array |I2-T^. Although the theoreti- 
cal investigation of this quantum phase transition has triggered 
enormous interest, the experimental realization of a large ar- 
ray of identical cavities is very ambitious. In this letter we 
show that even in the smallest possible coupled-cavity sys- 
tem of a Jaynes-Cummings dimer (JCD), one can find strong 
signatures of the on-site repulsive interaction among pho- 
tons. The two-coupled cavity system proposed here is sim- 
ple enough to be readily realizable with state-of-the-art circuit 
QED technology. 

We study a photon Josephson junction (PJJ) consisting of 
two tunnel-coupled microwave resonators each containing a 
single superconducting qubit (Fig.[Tl). We show that photons 
undergo a sharp non-equilibrium delocalization-localization 
transition from a regime where an initial photon population 
imbalance between the two resonators undergoes coherent 
oscillations (delocalized) between the two resonators to an- 
other regime where it becomes self-trapped (localized) as the 
photon-qubit interaction is increased. Similar self-trapping 
transitions were found in optical fibres Ii20|] . molecules 12 If . 
cold atom lEiUil and polariton BEC'sl25n. In all of these 
systems self-trapping is due to a Kerr/Bose-Hubbard like non- 



linearity and has been experimentally observed in the semi- 
classical regime with a large number of particles. The circuit 
QED implementation proposed in this paper has several ad- 
vantages with respect to these systems: (i) Self-trapping is due 
to a Jaynes-Cummings (JC) rather than a Kerr/Bose-Hubbard 
like nonlinearity. The JC interaction accurately describes the 
photon-qubit coupling in a microwave resonator fst] . In con- 
trast, a Kerr/Bose-Hubbard like nonlinearity is often a rather 
crude approximation of the experimental conditions 
(ii) The PJJ is a genuinely dissipative system. We show that 
dissipation and spontaneous emission of the qubit favor the 
localized regime, (iii) The PJJ may operate in the semiclas- 
sical (many photons) as well as quantum (few photons) limit 
since each resonator can initially be pumped with an almost 
arbitrary number of photons. 

In the following, we study in detail the classical versus quan- 
tum nature of this transition and present numerical as well as 
analytical results including the effects of dissipation and deco- 
herence. At the end of the paper, we outline a precise proposal 
on how to measure the localization transition of photons ex- 
perimentally. 

We describe the PJJ (Fig. [T]) by a two-site Jaynes- 
Cummings-Hubbard Hamiltonian (JCHM) 
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FIG. 1: Schematic diagram of the photonic Josephson junction (PJJ) 
proposed in this paper. Two transmission line microwave resonators 
(L,R) are coupled in series with a tunneling rate J, determined by 
the series capacitance of the resonators. Each resonator is strongly 
coupled to a superconducting qubit with a coupling rate g, providing 
a strong JC non-linearity. Photons can leave each resonator at a rate 
K, providing a mechanism for dissipation. Decoherence of the qubit 
is taken into account via the spontaneous emission rate 7. 
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FIG. 2: Photon imbalance z{t) and inversion of the qubits cr2(L.R)(t) 
obtained from a semiclassical approximation for the dissipative JCD 
with K — J = 0.05J. Shown are results in the Josephson regime 
with g — OAgc (left) and in the self-trapped regime with g = 2gc 
(right). The vertical dashed line points to a situation where all pho- 
tons have tunneled from the left to the right cavity and thus Rabi 
oscillations of the left qubit have slowed down while the right qubit 
oscillates rapidly. The inset in the right panel shows the inversion 
of the right qubit at short times in the self-trapped regime displaying 
strongly suppressed Rabi oscillations. Here, we have chosen the ini- 
tial condition ipL{0) ~ V^O corresponding to a coherent initial state 
with average photon number A'^ — 20. 



where /i^'-' denotes the local JC Hamiltonian hj'~^ — u!c a]ai + 
Lo^afa^ + g{a^ai + (J^al) for the left (L) or the right (R) 
cavity, ai (a|) and cr^ (c,^) are the photon creation (annihila- 
tion) and qubit raising (lowering) operators, respectively. The 
photon mode frequency is ojc, the qubit transition energy is uj^ 
and the photon-qubit coupling is given by g (we set h ^ 1). 
The photon-qubit interaction induces an anharmonicity in the 
spectrum of the JC Hamiltonian which leads to an effective 
on-site repulsion (anti-bunching) for photons. Throughout the 
paper we will assume zero detuning {uJx — (^c) for which this 
anharmonicity is strongest. Recently, dynamical aspects in a 
finite size JCHM have been investigated, but restricted to the 
one excitation (either photon or qubit) subspace (i.e., photon- 
photon interaction effects are irrelevant) and neglecting dissi- 
pation iBQ. 

In the JCD discussed here, photon dissipation and qubit de- 
cay are both taken into account by a Lindblad master equation 
for the system's density matrix p 



J2 i^C[a.]+^C[ar]), 
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where the Liouvillian of an operator O is defined as C[0] = 
(20/90t - Ot(9p _ pO^O) /2. Here, k and 7 denote cavity 
decay and spontaneous photon emission rate by the qubit, re- 
spectively, and we neglected pure qubit dephasing as it can be 
experimentally suppressed with respect to the other two decay 
channels, e.g. in a transmon qubit 13211 . Both decay processes 
in (O lead to a decrease of the total number of photons as a 
function of time. 
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FIG. 3: Photon imbalance z{t) obtained from a full numerical so- 
lution of the master equation l|2j (including quantum fluctuations) 
with K = 7 = (black thin curves) and n — j — 0.05J (red 
thick curves) for initially N{0) — 20 photons (and qubits initially in 
their groundstate (yz{'L,R){0) ~ —1/2). Figs, (a)-(f) show results for 
g = 0.1, 0.4, 0.6, 0.8, 0.9, 2gc. The inset shows the time-averaged 
imbalance (z) (averaged over the time interval t G [0, 100/ J]) for 
K = 7 = 0) (dashed line) as a function of the photon-qubit cou- 
pling g normalized with the semiclassical critical value gc in Eq. l[4j. 
In comparison, the semiclassical transition at this scale is essentially 
abrupt (dotted line). 

The central quantity of interest is the photon population im- 
balancez(t) = {7i-!^{t) — nY{{t))/N{t) withri^ — Tr al^aj/S and 
the total photon number = 71l + ^r. Throughout the pa- 
per, we consider the experimentally most relevant initial con- 
dition where the left cavity is pumped with N photons, the 
right cavity is empty and both qubits reside in their respec- 
tive ground-state att = 0. For large photon numbers, we can 
resort to a semiclassical approximation in which correlation 
functions in Eq. (|2]l are decoupled by simple factorization, 
e.g., (aV~) « (a^)(cr~), yielding eight coupled equations 
of motion for the expectation values of the photon and qubit 
operators which can be solved for an arbitrary number of pho- 
tons 113 311 . In the following, we present results of the semi- 
classical approximation (Fig. |2]l as well as the full quantum 
solution of Eq. (|2]i (Fig.[3]l- 

We first discuss results of the semiclassical approximation 
without dissipation (k = 7 = 0). In this case one can further 
reduce the number of coupled equations using energy conser- 
vation and the specific initial condition chosen above, yielding 
only four equations of motions 



Im(^R) 



-2gRe(VL) 
-2gIm(V'R,) 
.9siii(0L)/2- Jlm(V'R) 
.gsin(0R)/2 + JRe(V'L): 



(3) 



where V'(l.r) — (^(l.r))' Im(V'L) ~ R-e(?/'R) ~ and the an- 
gles 9i describe the qubit states (ctl) ~ — (sin6'L, 0, cos 0l)/2 
and ((Tr) = — (0, sin^R, cos6'r)/2. At zero interaction 
{g = 0) Eqs. ^ are exactly solvable. In this limit the 
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photon imbalance undergoes coherent harmonic oscillations 
z{t) = cosujjt with frequency ujj ~ 2 J. As g becomes non- 
zero this frequency decreases and the oscillations become an- 
harmonic. At a critical value of the coupling constant 



9c 



2.8VNJ 
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the period of oscillation diverges (critical slowing down) and 
an abrupt transition occurs to a localized regime, where the 
initial photon imbalance stays almost completely in the left 
cavity, i.e., z{t) w 1. Solutions of Eq. ([3]) near the transition 
(g ~ (7c) are shown in Fig.|2] An important and useful result 
is that the localization transition of photons can be observed 
in the population inversion of the two qubits, which depend 
on the number of photons in each cavity. In the delocalized 
regime (g < gc, see Fig. |2]left), when photons are tunneling, 
e.g., from the left into the right cavity, Rabi oscillations of 
the left qubit slow down considerably while those of the right 
qubit speed up. After half a tunneling period the scenario is 
reversed. On the other hand, in the localized regime (g > gc, 
see Fig. |2] right) the left qubit displays fast, complete Rabi 
oscillations while the right one displays slow, small amplitude 
oscillations (right column in Fig. |2]i. In other words, spatial 
localization of photons in the left resonator induces an energy 
state localization of the right qubit (notice that deep in the 
localized regime {g ^ gc) the right qubit remains very close 
to the ground-state at all times, i.e., (T^(l.r) {t) ~ Cz(l,r) (0)). 
This suggests that the localization transition of photons in a 
PJJ can be observed experimentally by a local readout of the 
qubit states. It also shows that qubit-qubit correlations are 
largely suppressed in the localized regime. 

Taking into account dissipation and spontaneous emission 
within the semiclassical approximation (i.e., solving all eight 
coupled equations of motion) leads to mainly three effects: 
(i) The oscillations of the qubit inversion are damped due to 
spontaneous emission, (ii) The frequency of these oscillations 
slow down with time as photons are leaving the cavity, (iii) 
Dissipation and spontaneous emission make the interaction 
terms of the Hamiltonian effectively time dependent. If the 
parameters of the system are properly chosen this can lead 
to a transition from a delocalized to a localized phase at fi- 
nite times. Thus, dissipation and spontaneous emission, sta- 
bilize the localized regime, i.e., we find self-trapped solutions 
for significantly smaller values of the photon-qubit coupling 
g < gc- In Fig.|2]we have chosen the coupling strengths such 
that we always stay within one phase; a careful analysis of this 
switching behavior will be given elsewhere. 

Full numerical solutions (including quantum fluctuations) 
of the master equation ^ using a Fock state basis for 20 
photons are shown in Fig. [3] Most importantly, we observe 
that the localization transition survives. However, specific 
quantum correction show up: (i) Deep inside the delocal- 
ized regime (see Fig. lUa)) the photon imbalance displays 
beatings of the coherent oscillations, a quantum feature that 
is absent in semiclassical solutions. Dissipation and spon- 
taneous emission strongly damp the large amplitude oscil- 
lations and suppress these beats; (ii) the localization tran- 
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FIG. 4: Rescaled photon imbalance z{t) = (ni(t) - n2(t))/iV(0) 
deep in the self-trapped regime for A'^(O) = 5 photons with g = S^c 
and n = 'y = Q. Shown are small amplitude Rabi oscillations on 
short time scales (Fig.|4ta)) and large amplitude ultra-long tunneling 
(Fig-Etb)). Results of a full numerical solution of the quantum mas- 
ter equation ([2J (grey curve in Figs.|4la) and|4lb)) are compared with 
strong-coupling degenerate perturbation theory (sdPT) (black curve) 
based on the effective level scheme shown in in Fig.|4jc). Note that 
a polariton eigenstate \M(j) is a mixed photon {M, M — 1) - qubit 
(g, e) state, i.e, |A/cr) = (|M , p) +cr|(M- 1) , e))/V2 (the zero po- 
lariton state is a special case with |0) = |0— ) = |0 , p)). In Fig.|4lc), 
Mcr, = {\Mo)L\Kfj)R, \Kfj)L\Mcj)R} denotes the pair of de- 
generate polariton eigenstates of the Hamiltonian (HJl at J = with 
{M, K) lower/upper (cr, /i = ±) polaritons in left (L) and right (R) 
cavity, respectively. Their degeneracy is lifted due to tunneling J, 
which induces the splitting A leading to ultra-long tunneling with 
period T = 2-k / /S. (see inset in Fig.|4tb)). 

sition is shifted to smaller g values and smoothened (see 
inset Fig. [3]); (iii) in the localized regime with dissipation 
and spontaneous emission (see Figs. [3le,f)) the imbalance 
approaches zero asymptotically at long times displaying no 
single-particle Rabi oscillations (i.e., when all photons left 
both cavities), (iv) the deeply localized regime displays rich 
multi-scale time dynamics, which can be explained using the 
effective level scheme shown in Fig. |4{c). For short times, 
Rabi oscillations of the qubit induce small amplitude oscilla- 
tions (on the order ~ 1 /N) of the rescaled photon imbalance 
z{t) = {ni{t) - n2{t))/N{0), see Fig.gta). The frequency 
of the Rabi oscillations (due to exchange of one photon be- 
tween local qubit and cavity) is given by the large splitting 
between lower and upper polariton states which yields for 
J = and large photon numbers ljb. = 2g\fN + ©(l/iV). 
For ultra-long times the localization of photons is unstable 
and almost complete oscillations of the imbalance set in, see 
Fig.Sfb). Neglecting the perturbative corrections to the eigen- 
vectors of the Hamiltonian ([T]i, we find for the rescaled imbal- 
ance z = cos(At) [l - (1/7V) sin2(w_Rt)] +0[J'^/g'^) which 
quantitatively reproduces the numerical results in Figs. 4(a) 
and 4(b) (red curve). The period of ultra-long tunneling is 
set by the splitting A, which we have calculated in leading 
order from A^-th order degenerate strong-coupling perturba- 
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tion theory, yielding A = CNJ{J/g)^~^ with a constant 
cpf that depends on the number of photons N. The inset in 
Fig. 4(b) shows the scaHng of the corresponding time period 
T = 27r/ A of the large amplitude oscillations as a function 
of N. The period T increases very fast (quasi-exponentially) 
with increasing photon number (see inset in Fig.|4|b)). Thus, 
true localization disappears in a small system, but turns expo- 
nentially 'good' with increasing system size (photon number). 
Already for five photons the ultra-long tunneling regime is 
hardly accessible experimentally. Ultra-long tunneling times 
were also predicted for the BHM 13111 . We should note, how- 
ever, that under experimental conditions any asymmetry be- 
tween the two cavities (e.g., due to small differences in de- 
tunings or coupling constants) will lift this degeneracy and 
set the effective time period of the large amplitude oscilla- 
tions. A corollary of this statement is that as long as such an 
asymmetry is much smaller than the frequency uij = 2 J the 
localization transition should be observable. 

The physics of a PJJ should be readily observable using 
a circuit QED implementation with realistic device parame- 
ters. A possible device consists of two series-coupled trans- 
mission line resonantors, each containing a single supercon- 
ducting qubit (Fig. [Til. A broad parameter space is avail- 
able through changes in lithographic patterning. In particular, 
qubit-cavity coupling g can range from 1 MHz to 300 MHz, 
while the cavity-cavity coupling J and the cavity dissipation 
rate k can be tuned independently in a range 50 kHz to 50 
MHz. The spontaneous emission rate 7 is typically 50KHz 
to 500kHz. An experimental observation of the localization 
transition proceeds in three parts. One cavity is populated 
with photons (initialization), the evolution proceeds for a fixed 
duration of time (evolution), and the photon occupancy of 
each cavity is finally measured (read-out). This entire process 
would be repeated for varying evolution times, thus allowing 
full reconstruction of the population imbalance z{t). Initial- 
ization can be accomplished using three different methods: (i) 
In the simplest method, the cavity is populated with a coherent 
photon state while the qubit is far off resonance {lUx <C uJc)', 
the qubit is then quickly brought into resonance (ujx = oJc) 
for the evolution 0, Etl- This scenario is best described by 
the results of the semiclassical approximation in Fig. |2] (ii) 
The cavity-qubit system can also be populated directly with a 
few-polariton state, i.e., an eigenstate of the JC Hamiltonian, 
using a properly timed 7r-pulse |6(| . (iii) Finally, a A^-photon 
Fock state can be constructed sequentially by successively ex- 
citing the qubit very quickly and bringing it into resonance 
(ujx — UJc). The multi-photon/polariton transitions in the lat- 
ter two scenarios are resolvable up to 5 — 10 excitations and 
are faithfully realized by the initial conditions chosen in Fig. [3] 
(full quantum calculation). After evolution for a given time, 
the qubits coupled to each cavity will be used to measure the 
photon occupation of that cavity. When strongly coupled, the 
qubit frequency is shifted depending on the number of photons 
in the cavity yl] ; by interrogating these different frequencies, 
a quantum non-demolition experiment of the photon number 



can be performed 10]. 

In this work we have shown that two tunnel-coupled mi- 
crowave resonators each containing a Jaynes-Cummings type 
nonlinearity undergo simultaneously sharp localization tran- 
sitions of photons (spatial) and qubits (energetic). Our results 
suggest many directions of further theoretical investigations 
including effects of detuning, quantum-classical crossover 
and interplay of localization and entanglement. 
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